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ABSTRACT 


This paper presents an analytical method of modelling eddy currents inside axial bearings. 
The problem is solved by dividing an axial bearing into elementary geometric forms, solving the 
Maxwell equations for these simplified geometries, defining boundary conditions and combining 
the geometries. The final result is an analytical solution for the flux, from which the impedance 
and the force of an axial bearing can be derived. Several impedance measurements have shown 
that the analytical solution can fit the measured data with a precision of approximately 5%. 


INTRODUCTION 


Modelling magnetic bearings is necessary to achieve reasonable results for magnetic force 
and bearing losses in order to design bearings and amplifiers. Furthermore, controller design 
is based on a model of the plant including magnetic bearings. Magnetic bearings are usu- 
ally modelled as equivalent electrical circuits. With simple models neglecting eddy currents, 
hysteresis, saturation and material nonlinearities, it is possible to achieve approximations for 
the magnetic force which are sufficient for most applications. Nevertheless, better modelling 
can improve bearing design and system behaviour. For some applications such as self sensing 

bearings a more precise model is necessary. . „ 

The absolute permeability //. describes the relationship between the magnetic flux density if 
and the magnetic field strength H (see equation (5)). The relative permeability // r is material 
dependent, the permeability of vacuum ji 0 is constant. For simple models p, can be assumed to 
be constant, but for real materials p. r depends on the magnetic field strength (i.e. // T - 
With increasing magnetic field strength the gradient /t decreases (to a limit value //. 0 ). There- 
fore the magnetic flux density is nearly constant for high field strengths ( saturation ). While 
increasing and decreasing H, the flux density B has different values for forward or backward 
loops due to hysteresis of the material. There are complex mathematical models such as the 
Preisach model [May91] to describe hysteresis, but these do not have an analytical solution. 
Changing magnetic fields inside conducting materials cause currents to flow. These eddy cur- 
rents have a large influence on the behaviour of magnetic bearings. The main goal of this paper 
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is to achieve an analytical description of a bearing model including eddy currents. Nonlinear 
material behaviour is considered in numerical calculations. Hysteresis is neglected in this paper. 
For elementary geometric forms, analytical solutions for eddy currents can be found. These 
solutions are sufficient to describe the behaviour of axial bearings. The elementary geometric 
forms are a semi-infinite plate, a rotationally symmetric plate and a semi-infinite cylinder. 


ELEMENTARY GEOMETRIC FORMS 


Electromagnetic fields and, therefore, eddy currents also can be described using Maxwell’s 
equations ([Kiip90], [Jac83]). In order to achieve analytical solutions it is necessary to divide 
an axial bearing into elementary geometric forms. Contrary to numerical results, analytical 
solutions can be more useful for bearing and controller design. Finite Element analysis can 
solve two or three-dimensional problems and more complex geometric forms. These numerical 
calculations are therefore used to verify analytical results and simplifications. 


(1) Semi-infinite cylinder 

(2) Rotational symmetric plate 

(3) Air gap 

(4) Coil with N turns 


Figure 1: Cross section of an axial bearing. The bearing is divided into elementary forms. 



Maxwell’s Equations 


Considering that the frequencies of the currents and fields are sufficiently low so that we can 
neglect the displacement current dD/dt, Maxwell’s equations can be written as 


curlJS = 


OB 

dt 


dB dH 

d H~dt 


curlif = J 


( 1 ) 

( 2 ) 


divJB = 0 


( 3 ) 


Here E is the electric field strength and J the current density. The generalized Ohm’s law, 
with the electric conductivity a is given by: 


J = oE 


( 4 ) 


The material equation is: 


B = fi 0 fi r H = flH 


( 5 ) 




Combining equations (1), (2), (4) and (5) we obtain: 


curl(curli/) = —fia 


OH 

1 H 


(G) 


A sinusoidal variation of the fields enables the transformation of equation (6) into the Fourier 
space, d/dt can be replaced by j u> giving: 

curl (curl//) = = —a 2 H (7) 

where: 

a =\/Mtv = ( 8 ) 

Maxwell’s equations (2) and (3) can be integrated using the theorems of Stokes and Gauss 
(©: magnetomotive force , <F: magnetic flux) to give: 

0 = J Hds = Y, I + j JdA (9) 

dA /I 

$ = J BdA (10) 

dv 

In the following sections, solutions for the elementary geometric forms semi-infinite plate , 
rotational symmetric plate and semi-infinite cylinder are given. 
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Figure 2: Cross section of 
(a) a semi-infinite plate, 



(b) a rotational plate, 


Hz(r) 

* J'P(r) 



(c) a semi-infinite cylinder 


Semi-infinite Plate 


A very long plate with thickness d and a magnetic field only in the z-dircction has an eddy 
current flow in the x-direction. Equation (7) leads to the following differential equation: 


( 0 

curl 0 

V H z 




(11) 
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a//z \ 
dy ' 

curl | 0 

0 ) 


0 

0 

-a 2 IL 


with the general solution: 


0 
0 

dy 2 

d 2 H z 2 

^ 


H z {y) = c x e ay + c 2 e- ny 


■J T {y) = -rr 1 = c x ac ay - c 2 ne oy 


dy 


( 12 ) 

(13) 

(14) 

(15) 


Rotational Symmetric Plate 


Hr 


0 \ 


curl | 0 | = | 

0 ) 


0 


= J, 


curl 


/ ° \ 

dilr 

dz 

0 / 


d 2 ll r 
' dz 2 

0 


-O' 


2 H r 


0 

0 


11 ( r dH A 

r dr V dz ) 

Equation (17) represents the following two differential equations: 


d 2 H r 2 


= a l H, 


dz 2 

1 dH r d 2 H r 


= 0 


(16) 


(17) 


(18) 


(19) 


r dz dzdr 

The solution for the differential equations (18) and (19) can be found by separation of the 
variables: 

H r (r,z) = 7 Z(r)Z(z) (20) 


Equation (18) can be transformed to: 


d 2 z 

dz 2 


= o?Z 


with the general solution: 


Z(z ) = z x e nz + z 2 e 
Equation (19) can be transformed to: 


1 dU n 
-77 + — — — 0 
r dr 


( 21 ) 

( 22 ) 

(23) 
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with the general solution: 


(24) 


n{r) - n- 

7 * 


Equation (22) and (24) lead to: 


1 


H r (r,z) = -( Cl e az +c 2 e~ nz ) 


(25) 


Without regard for r, equation (25) is similar to equation (14) (the case of the semi-infinite 
plate). 

Semi-infinite Cylinder 


curl 


l 0 

curl I 0 

V H z 


( 0 

dH z 

dr 
0 


/ 0 

an, 

dr 

o 


= Ju, 


0 

0 


\ -l £ (r^) 

\ r dr V dr J 


0 

0 

-q 2 H, 


(26) 


(27) 


d 2 H z 1 0H Z 


+ 


- a 2 H z = 0 


dr 2 r dr 

The substitution f = ar transforms the differential equation (28) into Bessel’s equation: 


(28) 




dr 2 


+ f- 


dr 


with the general solution: 


H z (r) = ci/ 0 (r) + c 2 I<o(r) 

or written without substitution: 

H z (r) - cJo(ar) + c 2 K 0 (ar) 


(29) 

(30) 

(31) 


Here, I 0 is the modified Bessel function of first type and zeroth order and I < 0 is the modified 
Bessel function of second type and zeroth order. More detailed information can be found in 
[AS65]. For calculations with Bessel functions, especially with complex arguments, a power 
series expansion is necessary. A solution without Bessel functions is therefore desirable. 

It can be clearly seen that equation (31) leads to equation (14) when » d (/?,: inner radius 
of the cylinder). A numerical comparison of the results for equation (31) and equation (14) 
shows that the difference between these two equations is only significant for unrealistically small 
values of Ri (Ri < d). 


An analysis of the cases of the semi-infinite cylinder and the rotational symmetric plate shows 
that these cases can be described with the equations of the semi-infinite plate. The problem of 
eddy currents in a semi-infinite plate is well known in the literature ([St.o74], [Kiip90]), because 
the choice of lamination thickness is important for high power transformers and electrical 
machines. 
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BOUNDARY CONDITIONS 


Solutions of differential equations require definitions of boundary conditions. Current field 
lines have to be closed. This is valid for eddy currents too. It is therefore assumed that eddy 
currents in a semi-infinite plate turn back at infinity. Eddy currents flow in one direction at 
one edge and in the other direction at the second edge (see figure 3 (a)). This solution is given 
in the literature ([Sto74], [Kiip90]). 

For a semi-infinite cylinder a second solution is possible. Contrary to a plate, a cylinder is 
geometrically closed and therefore eddy currents can flow in one direction (see figure 3 (b)). 
Calculations with the Finite Element program FEMAG ([Ins94]) produce this solution by de- 
fault. 



Case (a): A cylinder with an intersection. Case (b): A cylinder without an intersection. 
Eddy currents flow in two directions. Eddy currents flow in one direction. 

Figure 3: The two different boundary conditions. Cylinders with or without intersection. 


Cylinder with intersection 


The constants c x and c 2 of equations (14) and (15) can be found by defining two boundary 
conditions. The first boundary condition is given by 



T 

tola. 

II 

£ 

(32) 

The second boundary condition 

can be derived using Kirchhoff’s law. 

The sum of all the 

currents has to be 0. 

d 



2 

J J*(y)dy = 0 

d 

(33) 

Equation (33) can be written as: 

2 

I 1 




(34) 
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The solutions for H z , J T and the flux $ with the boundary condition of a cylinder with 
intersection are as follows (A is the area of the cross section of a cylinder): 


H z (y) = H 0 


cosli(ay) 

eo,sh(«|) 


(35) 


J x (y) = aH 0 


sinh(ary) 
cosh (a |) 


— j ^ d 

® = ~ / J* d Hz{y)dy = Ay, Ho 


tanh(al) 


*i 


(36) 


(37) 





0.1 Hz 
1 Hz 
10 Hz 


Figure 4: \H z (y)\ and \J x {y)\ for a cylinder with intersection. 

Figure 4 shows the decrease of the H - field in the middle of the core and the increase of eddy 
currents near the surface of the core with increasing frequency. The following model parameters 
have been used. 
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Cylinder without intersection 


In the case of a cylinder without intersection the start and end of the plate are fixed together 
and therefore eddy currents can have closed field lines and flow in one direction (the opposite 
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direction of the coil current). 1 The eddy current density at the outer surface of the core is zero. 


cj and c 2 in equations (14) and (15) can he found by defining two boundary conditions. The 
first boundary condition is given by: 

= H o (38) 

The second boundary condition is given by: 


j 4=° 


(39) 


The solutions for //*, 
are 


J x and $ with the boundary condition of a cylinder without intersection 


H,{y) = H 0 


cosh(a^ — ay) 
cosh(od) 


(40) 


J x {y) = aH 0 


sinh(a| — ay) 
cosh (ad) 


$ = AfiH 0 


tanh(od) 

ad 


(41) 

(42) 




y i m i xio - 3 


- / = 0.1 Hz 

- - / = 1 Hz 

- • / = 10 Hz 


Figure 5: \H z (y)\ and |Ji(y)| for a cylinder without intersection 

Figure 5 shows the decrease of the //-field near to the outer surface of the core and the 
increase of eddy currents near to the inner surface of the core with increasing frequency. The 
model parameters of the preceding section have been used. 


'The unidirectional eddy current flow reminds us of a transformer with a shorted secondary winding. With 
the transformer equation, a model for eddy currents would be very simple. A transformer however has two 
windings and one core. Currents flow inside the windings and the core is only used as a magnetic material. An 
axial bearing has one winding and the core has two functions (for conducting and as a magnetic material). The 
transformer equation therefore would lead to an unacceptable error. 
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Comparison of the two boundary conditions 


The generalized form of the flux equation can be written as: 

tanh(7) 

$ = $ 0 — 

7 


(43) 


with $ 0 = AhHq and with 7 = a\ . . . ad. 



f [Hz] 


Figure 6: |$(/)| 


cylinder without intersection 
cylinder with intersection 


A comparison of the flux equations (36) and (42) shows that 7 is two times larger for 
unidirectional than for bidirectional eddy currents. Figure 6 shows the frequency dependent 
flux for both cases with the same <J> 0 . 

The flux bandwidth of bidirectional eddy currents is 4 times higher (see figure 6). This does 
not correspond to measurements which only show a factor of 2. This difference is caused by 
the frequency dependence of 4>o- 

In the case of bidirectional eddy currents, 4> 0 is constant for all frequencies. The //-field only 
changes inside the core and the boundary conditions at the surface of the core remain the same. 
This leads to frequency independance of the global //-field. 

For unidirectional eddy currents the boundary conditions change at the surface of the core. 
This leads to a global change of the //-field and, therefore, to a change of <Fo- In order to 
achieve analytical solutions it is necessary to expand the problem of modelling eddy currents 
with the problem of modelling inductances. Calculation of inductances can be done for simple 
geometric setups, but for complex geometries such as axial bearings analytical models are not 
available. 

Several measurements have shown that the flux curves are similar for both eddy current cases. 
Therefore, the equations for bidirectional eddy currents can be used for unidirectional eddy 
currents too. It is clear that axial bearings with screw-holes or other gaps have uni- and 
bidirectional eddy currents and, therefore, the flux curve is between the two extreme eddy 
current cases. 
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ELECTRICAL MODEL AND FORCES 


The last boundary value H 0 can be found by closing the magnetic circuit. In order to build 
the complete magnetic circuit different materials have to be combined. Solutions for the flux 
are independent from the depth of the core (i.e. from the z-direction). Different materials lead 
to different solutions. It is clear that the flux changes in a transition region from one solution 
to the other. Finite Element analysis is used to calculate the flux distribution in a transition 
region between different materials. Figure 7 shows that the flux in the air gap is homogeneous 
even for high frequencies and that the transition region in the iron is small. Combining several 
materials is possible because the flux density distribution inside one material has little influence 
on the flux density distribution inside the next material. 



Figure 7: The flux inside the air gap is homogeneous even with a highly nonhomogeneous flux 
inside the unlaminated core (/ = lKHz). 


For calculations with eddy currents, it is useful to define the cut-off frequency f g of the flux 
decrease inside the core. In the literature f g is called the cut-off frequency of the lamination 
(see also [Kiip90]). 1 

fg = (44) 

Axial bearing materials normally have either a very low f g (unlaminated iron) or a very high 
f g (laminated iron, special magnetic core materials such as Corovac 2 * * and air ( f g<a ir —*■ °°))- 


Hdl = £ H„ 


= NI 


4> = / B m dA m = A m n 0 Hr,mH m tanh - — (46) 

J Tm 

2 Corovac is a high frequency core material with a « 5 ■ 10 5 S/in and /i r « 130. The /t r of Corovac is low 

compared to ferrite cores which are used for high frequency transformers. The machinability of Corovac is 

however much better. 
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Solving the system of equations (45) and (46) leads to the flux equation (for n different 
materials): N , 

* = — — ; ~ (47) 

y Im 

m=l Mr.rn A m tanh(7m) 

The modelling of axial bearings cannot be improved using more than two terms of the sum. 
The first term represents all materials with low f g and the second term includes all materials 
with a behaviour as air (/, is very high and tanh( 7 m )/ 7 m = 1). This simplification is based 
upon the fact that a connection of two materials with similar f g leads to: 


tanh( 7 i) tanh(72) 

a i — - + a 2 — - 

7i 72 


tanh( 7 3 ) 


The reduced flux equation is: 


VoNAl 

l hs. 7 

t^r tanh( 7 ) 


= HqNAI 


2x i | 1 7 

a fe tanh( 7 ) 


with: 


a fe Mr ‘ 


All areas are normalized to the area of the air gap A. x is the length of the air gap (including 
lengths of materials with very high f g ) and l fe is the length of the unlaminated core. 


Saturation 


In the previous calculations f.i r has been assumed to be constant. Saturation of the magnetic 
material leads to a limit value for the flux. The solution for the flux therefore gives unre- 
alistically high values at the surface. Even with simple models for saturation (e.g. Frohlich 
model [Sto74]), it is not possible to derive an analytical solution for the differential equations. 
Numerical calculations have shown that the behaviour of the H-field with saturation is very 
similar to the behaviour without saturation. It is therefore not necessary to include saturation 
into a model for eddy currents. 


Equivalent Electrical Circuit 


When x is assumed to be constant, the impedance of an eddy current affected coil can be 
written as: 

V A# 1 

Z =T= , =j»U- i (51) 


a fe tanh( 7 ) 


with: 


Lq = HqN 2 A- 
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Forces 


The magnetic force is proportional to the square of the flux. Therefore a reduction of the 
flux due to eddy currents will also reduce the force to: 


$ 2 _ KP 1 


(53) 


with: 


K = /t 0 N 2 A 


(54) 



Figure 8: Normalized force dynamics \F/F DC \ of an axial bearing (f g ~ 0.02Hz, a fe « 200) 


NONLINEAR PARAMETER IDENTIFICATION 


Parameter Identification 


It is clear that the impedance of a coil not only depends on the inductance affected by 
the eddy currents. An impedance model must also include the coil resistance R cu , the coil 
inductance L cu and the coil capacitance C cu (see figure 9). The last two parameters are necessary 
to model the frequency range above 10kHz. Thus, the total impedance is. 

I Rcu Lcu 

o— ► 

U 

▼ 

o — 

Figure 9: Complete equivalent electrical circuit 



432 




(55) 


Ztot — 


JojC c 


'll {Z + Rcu + jujLcu) 


Robustness of the identification can be improved by defining the constraints. The authors 
used the constraint algorithm from the MATLAB Optimization Toolbox. The parameter range 
described below is sufficient to approximate all measured axial bearings. Less stringent con- 
straints are possible, but computational time increases. 

All parameter ranges excluding R cu are written with exponential notation. The fitting al- 
gorithm therefore fits the exponent of the parameters. This method avoids the value 0 and 
guarantees convergence. 


The ranges of the fitted parameters are given below. 
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Figure 10: Example 1: Axial bearing with unlaminated core and target. 

— : calculated data, - - : measured data 

Identified parameters: R cu « C cu « 800pF, L cu ss 120 iiH, L 0 ~ 12 mH, Of e ~ 3, f g « 711 z 
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Figure 11: Example 2: Axial bearing with a Corovac core and an unlaminated target. 

— : calculated data, - - : measured data 

Identified parameters: R cu ~ 0.4Q, Cc ~ 150 pF,L cu ~ Lq ~ \AmH,df e w 150, f g ~ 

0.12 Hz 


The fitting error e is approximately 5% with the usual air gap and increases up to 10% when 
the air gap is small and the flux is saturated and is given by: 




Zfit 

Z-meaa 


2 


(56) 


where n is the number of measured values, Z fit the fitted impedance and Z meaa the measured 
impedance. 


Losses 


The spectral losses of an axial bearing are given by: 

P = real(t/conj(/)) = |f/| 2 real(conj )) (57) 

\Ztot' 

Figure 12 shows normalized loss based on impedance measurements of an axial bearing with 
a Corovac core. At low frequencies losses are very high due to R cu . With increasing frequency 
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losses decrease. It is therefore useful to choose the switching frequency of an amplifier as high 
as possible. As mentioned in section Comparison of the two boundary conditions, a radial 
intersection of an axial bearing can increase the bandwidth of the flux and leads to lower losses. 
It is only sufficient to intersect unlaminated parts of the core with low f g , i.e. in most cases the 
axial bearing target. An axial bearing without intersection suffers more than 40% higher loss 
than one with intersection at frequencies typically used for switching amplifiers (20 - 200 kHz). 




Figure 12: upper plot: normalized losses P n = P/\n \ 2 of an axial bearing with Corovac. core 

and steel target (— P x \ target intersected , - - P 2 \ target not intersected ) 

lower plot: relative difference between the two curves of the upper plot. P r — (P 2 - P\)/P\ 

Measurement 


Measurements of axial bearing impedances are made with a LCR-meter (HP 4284A) includ- 
ing a power current source (HP 42841 A). The LCR-meter has the possibility to measure with 
logarithmically sweeped frequencies from 20Hz up to 1MHz. A host computer controls the 
measurement via GPIB-interface. 

The power current source allows measurements of very high inductances and can superpose the 
measurement signal with a DC-bias current of up to 20A. 


CONCLUSIONS 


A model for eddy currents in axial bearings has been derived which can describe the be- 
haviour for a wide frequency range. From the analytical solution for the flux, the impedance, 
the magnetic force and the losses can be calculated. Measurements have shown that the model 
is very precise. 
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